{smcl}
{com}{sf}{ul off}{txt}{.-}
      name:  {res}<unnamed>
       {txt}log:  {res}C:\Users\jp18390\Dropbox\DISCRIMINATION PROJECT (U.S. FEDERAL AGENCIES)\Age Discrimination Project\PAR R&R\Statistics\Age_Discrimination.APPENDIX 3.ONLY LONG PANELS.06-02-2022.smcl
  {txt}log type:  {res}smcl
 {txt}opened on:  {res} 2 Jun 2022, 18:09:51
{txt}
{com}.    
.    
. 
. use "C:\Users\jp18390\Dropbox\DISCRIMINATION PROJECT (U.S. FEDERAL AGENCIES)\Age Discrimination Project\PAR R&R\Statistics\Age_Discrimination.Post-Estimation.06-02-2022.dta", replace 
{txt}
{com}. 
. 
. 
. 
. xtset a_id year, yearly
{res}{txt}{col 8}panel variable:  {res}a_id (unbalanced)
{txt}{col 9}time variable:  {res}{col 25}year, 2010 to 2019, but with gaps
{txt}{col 17}delta:  {res}1 year
{txt}
{com}. 
.   
. 
. **** 2022 PAR Data Replication [06/02/2022] "UNDERSTANDING ORGANIZATIONAL SUSCEPTIBILITY TO AGE DISCIMINRATION WITHIN THE U.S. FEDERAL WORKFORCE" [KRAUSE & PARK] ****
. 
. 
. 
.  
. **** APPENDIX SECTION 3: PANEL SENSITIVITY ANALYSES -- INCLUSION OF ONLY AGENCIES WITH AT LEAST SEVEN ANNUAL OBSERVATIONS [T-BAR = 6.9]  **** 
. 
.    
.    
. *****************************************************************************************************************************************************************************************
. *****************************************************************************************************************************************************************************************
. *****************************************************************************************************************************************************************************************
. *****************************************************************************************************************************************************************************************
.    
. *** GENERATE IDENTIFER FOR PANELS OF VARYING SIZES: FULL PANELS [T=10] = 2; WEAKLY UNBALANCED PANELS [9 =< T >= 7: ABOVE THE MEAN ANUAL OBSERVATIONS = 6.9]; SHORTER PANELS [T < 7] ***
. 
. *** 38 agencies contain full panels of observations [T= 10] ***
. gen suffpaneldum = 2 if a_id==1 | a_id==8 | a_id==9 | a_id==10 |  a_id== 35 | a_id==38 | a_id==40 | a_id==44 |  a_id== 50 | a_id==51 | a_id==53 | a_id==54 | a_id==55 | a_id==56 | a_id==57 | a_id==58 |  a_id== 59 | a_id==61 | a_id==62 | a_id==65 |  a_id==66 | a_id==75 | a_id==76 | a_id==77 | a_id==79 | a_id==80 | a_id==81 | a_id==83 | a_id==91 | a_id==92 | a_id==96 | a_id==102 |  a_id==112 | a_id==115 | a_id==123 | a_id==124 |  a_id== 127 | a_id==130
{txt}(517 missing values generated)

{com}. 
. *** 43 agencies contain 7-9 observations per panel [9 <= T >= 7] ***
. 
. replace suffpaneldum = 1 if a_id==6 | a_id==7 | a_id==12 | a_id==13 | a_id== 14 | a_id==31 | a_id==32 | a_id==33 |  a_id==34 | a_id==37 | a_id==39 | a_id==41 | a_id==42 | a_id==43 | a_id==48 | a_id==49 | a_id==52 | a_id==67 | a_id==68 | a_id==69 | a_id==70 | a_id==71 | a_id==72 | a_id==73 | a_id==74 | a_id==82 | a_id==85 | a_id==87 | a_id==94 | a_id==95 | a_id==97 | a_id==100 | a_id==101 | a_id==103 | a_id==105 | a_id==108 | a_id== 111 | a_id==116 | a_id==117 | a_id==120 | a_id==121 | a_id==125 | a_id==129
{txt}(344 real changes made)

{com}. 
. *** 49 agencies contain Less than 7 observations per panel [T < 7] ***
. replace suffpaneldum = 0 if suffpaneldum==.
{txt}(173 real changes made)

{com}. 
. tab suffpaneldum 

{txt}suffpaneldu {c |}
          m {c |}      Freq.     Percent        Cum.
{hline 12}{c +}{hline 35}
          0 {c |}{res}        173       19.29       19.29
{txt}          1 {c |}{res}        344       38.35       57.64
{txt}          2 {c |}{res}        380       42.36      100.00
{txt}{hline 12}{c +}{hline 35}
      Total {c |}{res}        897      100.00
{txt}
{com}. 
. 
. 
. 
. 
. 
. 
. *** TABLE A3.2: REGRESSION MODEL TABLE PREDICTING VARIATIONS IN AGE DISCRIMINATION FORMAL COMPLAINTS IN U.S. FEDERAL AGENCIES ***
. 
. 
. 
. 
. 
. ** MODEL A3.1: DISAGGREGATE SUPERVISOR/SUBORDINATE RATIO MEASURES [ratio_40over_suplb; ratio_40over_nonsuplb]: ONLY RANDOM INTERCEPT MODEL SPECIFICATION WITH BE & WE ESTIMATES FOR ALL COVARIATES [SANS YEAR UNIT EFFECTS & LN(TOTAL WORKFORCE)] --- REDUCED MODEL: ONLY CONTROL COVARIATES ARE LNTEF & YEAR UNIT EFFECTS ***
. 
. xthybrid age_discrimination  ratio_40over_suplb  ratio_40over_nonsuplb   lntwf ///
> yr2 yr3 yr4 yr5 yr6 yr7 yr8 yr9 yr10 ratio_55over_totallb if suffpaneldum!=0 , clusterid(a_id) vce(cluster a_id)  family(nbinomial) link(log) full ///
> use(ratio_40over_suplb  ratio_40over_nonsuplb ratio_55over_totallb) 

{res}
{txt}{hline}
{p 0 8}Model {hi:model}{p_end}
{hline}

{txt}Mixed-effects GLM{col 49}{txt}Number of obs{col 67}={res}{col 69}       724
{txt}Family: {col 15}{res}negative binomial
{txt}Link: {col 29}{res}log
{txt}Overdispersion: {col 28}{res}mean
{txt}Group variable: {col 28}{res}a_id{col 49}{txt}Number of groups{col 67}={res}{col 69}        81

{col 49}{txt}Obs per group:
{col 63}{txt}min{col 67}={res}{col 69}         7
{col 63}{txt}avg{col 67}={res}{col 69}       8.9
{col 63}{txt}max{col 67}={res}{col 69}        10

{txt}Integration method: {col 21}{res}mvaghermite{col 49}{txt}Integration pts.{col 67}={col 78}{res}7

{col 49}{txt}Wald chi2({res}16{txt}){col 67}={res}{col 70}   761.27
{txt}Log pseudolikelihood = {res}-2348.7731{col 49}{txt}Prob > chi2{col 67}={res}{col 73}0.0000
{txt}{ralign 78:(Std. Err. adjusted for {res:81} clusters in a_id)}
{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 14}{c |}{col 26}    Robust
{col 1}age_discri~n{col 14}{c |}      Coef.{col 26}   Std. Err.{col 38}      z{col 46}   P>|z|{col 54}     [95% Con{col 67}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 6}R__yr2 {c |}{col 14}{res}{space 2}-.0212766{col 26}{space 2} .0508754{col 37}{space 1}   -0.42{col 46}{space 3}0.676{col 54}{space 4}-.1209905{col 67}{space 3} .0784373
{txt}{space 6}R__yr3 {c |}{col 14}{res}{space 2}-.0663139{col 26}{space 2} .0587355{col 37}{space 1}   -1.13{col 46}{space 3}0.259{col 54}{space 4}-.1814333{col 67}{space 3} .0488055
{txt}{space 6}R__yr4 {c |}{col 14}{res}{space 2}-.1152551{col 26}{space 2}  .068373{col 37}{space 1}   -1.69{col 46}{space 3}0.092{col 54}{space 4}-.2492637{col 67}{space 3} .0187535
{txt}{space 6}R__yr5 {c |}{col 14}{res}{space 2}-.1237208{col 26}{space 2} .0810699{col 37}{space 1}   -1.53{col 46}{space 3}0.127{col 54}{space 4}-.2826149{col 67}{space 3} .0351733
{txt}{space 6}R__yr6 {c |}{col 14}{res}{space 2} -.015667{col 26}{space 2} .0746022{col 37}{space 1}   -0.21{col 46}{space 3}0.834{col 54}{space 4}-.1618847{col 67}{space 3} .1305506
{txt}{space 6}R__yr7 {c |}{col 14}{res}{space 2}-.0070046{col 26}{space 2} .0762397{col 37}{space 1}   -0.09{col 46}{space 3}0.927{col 54}{space 4}-.1564316{col 67}{space 3} .1424224
{txt}{space 6}R__yr8 {c |}{col 14}{res}{space 2}-.0994051{col 26}{space 2} .0770461{col 37}{space 1}   -1.29{col 46}{space 3}0.197{col 54}{space 4}-.2504127{col 67}{space 3} .0516026
{txt}{space 6}R__yr9 {c |}{col 14}{res}{space 2}-.0596529{col 26}{space 2} .0790809{col 37}{space 1}   -0.75{col 46}{space 3}0.451{col 54}{space 4}-.2146486{col 67}{space 3} .0953428
{txt}{space 5}R__yr10 {c |}{col 14}{res}{space 2}-.0816361{col 26}{space 2}  .091319{col 37}{space 1}   -0.89{col 46}{space 3}0.371{col 54}{space 4} -.260618{col 67}{space 3} .0973458
{txt}{space 4}R__lntwf {c |}{col 14}{res}{space 2} .7902359{col 26}{space 2} .0369974{col 37}{space 1}   21.36{col 46}{space 3}0.000{col 54}{space 4} .7177223{col 67}{space 3} .8627495
{txt}W__ra~_suplb {c |}{col 14}{res}{space 2} .0768014{col 26}{space 2} 1.159586{col 37}{space 1}    0.07{col 46}{space 3}0.947{col 54}{space 4}-2.195945{col 67}{space 3} 2.349547
{txt}W__ra~nsuplb {c |}{col 14}{res}{space 2} .3243907{col 26}{space 2} .9414957{col 37}{space 1}    0.34{col 46}{space 3}0.730{col 54}{space 4}-1.520907{col 67}{space 3} 2.169688
{txt}W__ratio_5~b {c |}{col 14}{res}{space 2} 1.896918{col 26}{space 2} 1.022162{col 37}{space 1}    1.86{col 46}{space 3}0.063{col 54}{space 4} -.106483{col 67}{space 3}  3.90032
{txt}B__ra~_suplb {c |}{col 14}{res}{space 2}-5.047685{col 26}{space 2} .9572618{col 37}{space 1}   -5.27{col 46}{space 3}0.000{col 54}{space 4}-6.923884{col 67}{space 3}-3.171487
{txt}B__ra~nsuplb {c |}{col 14}{res}{space 2} 2.777659{col 26}{space 2} 1.097347{col 37}{space 1}    2.53{col 46}{space 3}0.011{col 54}{space 4} .6268983{col 67}{space 3} 4.928421
{txt}B__ratio_5~b {c |}{col 14}{res}{space 2}-.0668453{col 26}{space 2} 1.053153{col 37}{space 1}   -0.06{col 46}{space 3}0.949{col 54}{space 4}-2.130988{col 67}{space 3} 1.997297
{txt}{space 7}_cons {c |}{col 14}{res}{space 2}-1.740413{col 26}{space 2} .8427916{col 37}{space 1}   -2.07{col 46}{space 3}0.039{col 54}{space 4}-3.392254{col 67}{space 3}-.0885718
{txt}{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 4}/lnalpha {c |}{col 14}{res}{space 2}-3.299953{col 26}{space 2} .1805961{col 54}{space 4}-3.653914{col 67}{space 3}-2.945991
{txt}{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 1}{res}a_id        {col 14}{txt}{c |}
   var(_cons){c |}{col 14}{res}{space 2} .1633559{col 26}{space 2} .0294498{col 54}{space 4} .1147308{col 67}{space 3} .2325893
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}

{com}. *
. estat ic

Akaike's information criterion and Bayesian information criterion

{txt}{hline 13}{c TT}{hline 63}
       Model {c |}        Obs  ll(null)  ll(model)      df         AIC        BIC
{hline 13}{c +}{hline 63}
{ralign 12:{stata estimates replay model:model}}{col 14}{c |}{res}{col 16}       724{col 27}        .{col 38}-2348.773{col 49}    19{col 58} 4735.546{col 69} 4822.657
{txt}{hline 13}{c BT}{hline 63}
{p 15 21 2}
Note: N=Obs used in calculating BIC; see {helpb bic_note:[R] BIC note}.
{p_end}

{com}. 
. 
. *
. *
. *
. *
. *
. 
. ** MODEL A3.2: DISAGGREGATE SUPERVISOR/SUBORDINATE RATIO MEASURES [ratio_40over_suplb; ratio_40over_nonsuplb]: ONLY RANDOM INTERCEPT MODEL SPECIFICATION WITH BE & WE ESTIMATES FOR ALL COVARIATES [SANS YEAR UNIT EFFECTS & LN(TOTAL WORKFORCE)] --- FULL MODEL: ALL CONTROL COVARIATES ARE LNTEF & YEAR UNIT EFFECTS ***
. 
. xthybrid age_discrimination  ratio_40over_suplb  ratio_40over_nonsuplb   orgjustice_sem   nonprof40over_tr_lb   politicization_lb ratio_fsup_msup ratio_minsup_nonmsup lntwf ///
> yr2 yr3 yr4 yr5 yr6 yr7 yr8 yr9 yr10 ratio_55over_totallb  if suffpaneldum!=0, clusterid(a_id) vce(cluster a_id)  family(nbinomial) link(log) full ///
> use(ratio_40over_suplb  ratio_40over_nonsuplb    orgjustice_sem  nonprof40over_tr_lb   politicization_lb   ratio_fsup_msup ratio_minsup_nonmsup ratio_55over_totallb) 

{res}
{txt}{hline}
{p 0 8}Model {hi:model}{p_end}
{hline}

{txt}Mixed-effects GLM{col 49}{txt}Number of obs{col 67}={res}{col 69}       724
{txt}Family: {col 15}{res}negative binomial
{txt}Link: {col 29}{res}log
{txt}Overdispersion: {col 28}{res}mean
{txt}Group variable: {col 28}{res}a_id{col 49}{txt}Number of groups{col 67}={res}{col 69}        81

{col 49}{txt}Obs per group:
{col 63}{txt}min{col 67}={res}{col 69}         7
{col 63}{txt}avg{col 67}={res}{col 69}       8.9
{col 63}{txt}max{col 67}={res}{col 69}        10

{txt}Integration method: {col 21}{res}mvaghermite{col 49}{txt}Integration pts.{col 67}={col 78}{res}7

{col 49}{txt}Wald chi2({res}26{txt}){col 67}={res}{col 70}  1067.49
{txt}Log pseudolikelihood = {res} -2330.905{col 49}{txt}Prob > chi2{col 67}={res}{col 73}0.0000
{txt}{ralign 78:(Std. Err. adjusted for {res:81} clusters in a_id)}
{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 14}{c |}{col 26}    Robust
{col 1}age_discri~n{col 14}{c |}      Coef.{col 26}   Std. Err.{col 38}      z{col 46}   P>|z|{col 54}     [95% Con{col 67}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 6}R__yr2 {c |}{col 14}{res}{space 2}-.0115822{col 26}{space 2} .0495244{col 37}{space 1}   -0.23{col 46}{space 3}0.815{col 54}{space 4}-.1086483{col 67}{space 3} .0854838
{txt}{space 6}R__yr3 {c |}{col 14}{res}{space 2}-.0500762{col 26}{space 2}  .054992{col 37}{space 1}   -0.91{col 46}{space 3}0.363{col 54}{space 4}-.1578586{col 67}{space 3} .0577062
{txt}{space 6}R__yr4 {c |}{col 14}{res}{space 2} -.096633{col 26}{space 2} .0672403{col 37}{space 1}   -1.44{col 46}{space 3}0.151{col 54}{space 4}-.2284215{col 67}{space 3} .0351556
{txt}{space 6}R__yr5 {c |}{col 14}{res}{space 2}-.0988012{col 26}{space 2} .0808048{col 37}{space 1}   -1.22{col 46}{space 3}0.221{col 54}{space 4}-.2571756{col 67}{space 3} .0595732
{txt}{space 6}R__yr6 {c |}{col 14}{res}{space 2}  .011793{col 26}{space 2} .0759254{col 37}{space 1}    0.16{col 46}{space 3}0.877{col 54}{space 4} -.137018{col 67}{space 3}  .160604
{txt}{space 6}R__yr7 {c |}{col 14}{res}{space 2} .0148696{col 26}{space 2} .0775068{col 37}{space 1}    0.19{col 46}{space 3}0.848{col 54}{space 4}-.1370409{col 67}{space 3} .1667801
{txt}{space 6}R__yr8 {c |}{col 14}{res}{space 2}-.0865635{col 26}{space 2}  .079608{col 37}{space 1}   -1.09{col 46}{space 3}0.277{col 54}{space 4}-.2425923{col 67}{space 3} .0694652
{txt}{space 6}R__yr9 {c |}{col 14}{res}{space 2}-.0350252{col 26}{space 2}  .080141{col 37}{space 1}   -0.44{col 46}{space 3}0.662{col 54}{space 4}-.1920986{col 67}{space 3} .1220483
{txt}{space 5}R__yr10 {c |}{col 14}{res}{space 2}-.0519484{col 26}{space 2} .0920595{col 37}{space 1}   -0.56{col 46}{space 3}0.573{col 54}{space 4}-.2323817{col 67}{space 3} .1284849
{txt}{space 4}R__lntwf {c |}{col 14}{res}{space 2} .7943112{col 26}{space 2} .0367391{col 37}{space 1}   21.62{col 46}{space 3}0.000{col 54}{space 4} .7223039{col 67}{space 3} .8663186
{txt}W__ra~_suplb {c |}{col 14}{res}{space 2} -.008714{col 26}{space 2} 1.188766{col 37}{space 1}   -0.01{col 46}{space 3}0.994{col 54}{space 4}-2.338652{col 67}{space 3} 2.321224
{txt}W__ra~nsuplb {c |}{col 14}{res}{space 2} .9181165{col 26}{space 2} 1.160585{col 37}{space 1}    0.79{col 46}{space 3}0.429{col 54}{space 4}-1.356588{col 67}{space 3} 3.192821
{txt}W__orgjust~m {c |}{col 14}{res}{space 2}  .157321{col 26}{space 2}  .257217{col 37}{space 1}    0.61{col 46}{space 3}0.541{col 54}{space 4} -.346815{col 67}{space 3}  .661457
{txt}W__nonprof~b {c |}{col 14}{res}{space 2}-.7192953{col 26}{space 2} 1.437199{col 37}{space 1}   -0.50{col 46}{space 3}0.617{col 54}{space 4}-3.536154{col 67}{space 3} 2.097563
{txt}W__politic~b {c |}{col 14}{res}{space 2}-1.880088{col 26}{space 2} 1.670104{col 37}{space 1}   -1.13{col 46}{space 3}0.260{col 54}{space 4} -5.15343{col 67}{space 3} 1.393255
{txt}W__ratio_f~p {c |}{col 14}{res}{space 2} .2500767{col 26}{space 2} .2513558{col 37}{space 1}    0.99{col 46}{space 3}0.320{col 54}{space 4}-.2425716{col 67}{space 3} .7427251
{txt}W__ratio_m~p {c |}{col 14}{res}{space 2}-.3566136{col 26}{space 2} .1180265{col 37}{space 1}   -3.02{col 46}{space 3}0.003{col 54}{space 4}-.5879412{col 67}{space 3}-.1252859
{txt}W__ratio_5~b {c |}{col 14}{res}{space 2} 1.930509{col 26}{space 2} 1.029501{col 37}{space 1}    1.88{col 46}{space 3}0.061{col 54}{space 4}-.0872755{col 67}{space 3} 3.948293
{txt}B__ra~_suplb {c |}{col 14}{res}{space 2}    -3.59{col 26}{space 2} .9718747{col 37}{space 1}   -3.69{col 46}{space 3}0.000{col 54}{space 4} -5.49484{col 67}{space 3}-1.685161
{txt}B__ra~nsuplb {c |}{col 14}{res}{space 2} 2.272572{col 26}{space 2} .9630955{col 37}{space 1}    2.36{col 46}{space 3}0.018{col 54}{space 4} .3849394{col 67}{space 3} 4.160204
{txt}B__orgjust~m {c |}{col 14}{res}{space 2}-1.343862{col 26}{space 2} .3508769{col 37}{space 1}   -3.83{col 46}{space 3}0.000{col 54}{space 4}-2.031568{col 67}{space 3}-.6561558
{txt}B__nonprof~b {c |}{col 14}{res}{space 2} .1390077{col 26}{space 2} .2673495{col 37}{space 1}    0.52{col 46}{space 3}0.603{col 54}{space 4}-.3849878{col 67}{space 3} .6630031
{txt}B__politic~b {c |}{col 14}{res}{space 2} 1.703256{col 26}{space 2} .7990719{col 37}{space 1}    2.13{col 46}{space 3}0.033{col 54}{space 4} .1371039{col 67}{space 3} 3.269408
{txt}B__ratio_f~p {c |}{col 14}{res}{space 2} .1609125{col 26}{space 2} .0903123{col 37}{space 1}    1.78{col 46}{space 3}0.075{col 54}{space 4}-.0160963{col 67}{space 3} .3379213
{txt}B__ratio_m~p {c |}{col 14}{res}{space 2}-.0407889{col 26}{space 2} .0509883{col 37}{space 1}   -0.80{col 46}{space 3}0.424{col 54}{space 4}-.1407242{col 67}{space 3} .0591464
{txt}B__ratio_5~b {c |}{col 14}{res}{space 2}-.3797641{col 26}{space 2}  .829731{col 37}{space 1}   -0.46{col 46}{space 3}0.647{col 54}{space 4}-2.006007{col 67}{space 3} 1.246479
{txt}{space 7}_cons {c |}{col 14}{res}{space 2}-2.740205{col 26}{space 2} .8476294{col 37}{space 1}   -3.23{col 46}{space 3}0.001{col 54}{space 4}-4.401529{col 67}{space 3}-1.078882
{txt}{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 4}/lnalpha {c |}{col 14}{res}{space 2}-3.340742{col 26}{space 2} .1750734{col 54}{space 4}-3.683879{col 67}{space 3}-2.997604
{txt}{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 1}{res}a_id        {col 14}{txt}{c |}
   var(_cons){c |}{col 14}{res}{space 2} .1216228{col 26}{space 2} .0264472{col 54}{space 4} .0794179{col 67}{space 3} .1862564
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}

{com}. *
. estat ic

Akaike's information criterion and Bayesian information criterion

{txt}{hline 13}{c TT}{hline 63}
       Model {c |}        Obs  ll(null)  ll(model)      df         AIC        BIC
{hline 13}{c +}{hline 63}
{ralign 12:{stata estimates replay model:model}}{col 14}{c |}{res}{col 16}       724{col 27}        .{col 38}-2330.905{col 49}    29{col 58}  4719.81{col 69} 4852.769
{txt}{hline 13}{c BT}{hline 63}
{p 15 21 2}
Note: N=Obs used in calculating BIC; see {helpb bic_note:[R] BIC note}.
{p_end}

{com}. 
. 
. 
. 
. 
. 
. 
. ** MODEL A3.3: RATIO OF OVER40 SUPERVISORS TO OVER40 NON-SUPERVISORS [RATIO OF ratio_40over_suplb TO ratio_40over_nonsuplb]: ONLY RANDOM INTERCEPT MODEL SPECIFICATION WITH BE & WE ESTIMATES FOR ALL COVARIATES [SANS YEAR UNIT EFFECTS & LN(TOTAL WORKFORCE)] REDUCED MODEL: ONLY CONTROL COVARIATES ARE LNTEF & YEAR UNIT EFFECTS ***
. 
. xthybrid age_discrimination  ratio40suplb_nonsuplb   lntwf ///
> yr2 yr3 yr4 yr5 yr6 yr7 yr8 yr9 yr10 ratio_55over_totallb if suffpaneldum!=0 , clusterid(a_id) vce(cluster a_id)  family(nbinomial) link(log) full ///
> use(ratio40suplb_nonsuplb ratio_55over_totallb) 

{res}
{txt}{hline}
{p 0 8}Model {hi:model}{p_end}
{hline}

{txt}Mixed-effects GLM{col 49}{txt}Number of obs{col 67}={res}{col 69}       724
{txt}Family: {col 15}{res}negative binomial
{txt}Link: {col 29}{res}log
{txt}Overdispersion: {col 28}{res}mean
{txt}Group variable: {col 28}{res}a_id{col 49}{txt}Number of groups{col 67}={res}{col 69}        81

{col 49}{txt}Obs per group:
{col 63}{txt}min{col 67}={res}{col 69}         7
{col 63}{txt}avg{col 67}={res}{col 69}       8.9
{col 63}{txt}max{col 67}={res}{col 69}        10

{txt}Integration method: {col 21}{res}mvaghermite{col 49}{txt}Integration pts.{col 67}={col 78}{res}7

{col 49}{txt}Wald chi2({res}14{txt}){col 67}={res}{col 70}   637.77
{txt}Log pseudolikelihood = {res}-2352.6387{col 49}{txt}Prob > chi2{col 67}={res}{col 73}0.0000
{txt}{ralign 78:(Std. Err. adjusted for {res:81} clusters in a_id)}
{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 14}{c |}{col 26}    Robust
{col 1}age_discri~n{col 14}{c |}      Coef.{col 26}   Std. Err.{col 38}      z{col 46}   P>|z|{col 54}     [95% Con{col 67}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 6}R__yr2 {c |}{col 14}{res}{space 2}-.0248236{col 26}{space 2} .0516223{col 37}{space 1}   -0.48{col 46}{space 3}0.631{col 54}{space 4}-.1260014{col 67}{space 3} .0763543
{txt}{space 6}R__yr3 {c |}{col 14}{res}{space 2}-.0717714{col 26}{space 2} .0607511{col 37}{space 1}   -1.18{col 46}{space 3}0.237{col 54}{space 4}-.1908413{col 67}{space 3} .0472985
{txt}{space 6}R__yr4 {c |}{col 14}{res}{space 2}-.1238916{col 26}{space 2} .0704024{col 37}{space 1}   -1.76{col 46}{space 3}0.078{col 54}{space 4}-.2618777{col 67}{space 3} .0140946
{txt}{space 6}R__yr5 {c |}{col 14}{res}{space 2}-.1356444{col 26}{space 2} .0825602{col 37}{space 1}   -1.64{col 46}{space 3}0.100{col 54}{space 4}-.2974594{col 67}{space 3} .0261707
{txt}{space 6}R__yr6 {c |}{col 14}{res}{space 2}-.0307759{col 26}{space 2} .0776961{col 37}{space 1}   -0.40{col 46}{space 3}0.692{col 54}{space 4}-.1830575{col 67}{space 3} .1215056
{txt}{space 6}R__yr7 {c |}{col 14}{res}{space 2}-.0236774{col 26}{space 2} .0796292{col 37}{space 1}   -0.30{col 46}{space 3}0.766{col 54}{space 4}-.1797478{col 67}{space 3}  .132393
{txt}{space 6}R__yr8 {c |}{col 14}{res}{space 2}-.1166431{col 26}{space 2} .0796699{col 37}{space 1}   -1.46{col 46}{space 3}0.143{col 54}{space 4}-.2727933{col 67}{space 3}  .039507
{txt}{space 6}R__yr9 {c |}{col 14}{res}{space 2}-.0779929{col 26}{space 2} .0820497{col 37}{space 1}   -0.95{col 46}{space 3}0.342{col 54}{space 4}-.2388074{col 67}{space 3} .0828215
{txt}{space 5}R__yr10 {c |}{col 14}{res}{space 2}-.0991713{col 26}{space 2} .0937317{col 37}{space 1}   -1.06{col 46}{space 3}0.290{col 54}{space 4}-.2828821{col 67}{space 3} .0845395
{txt}{space 4}R__lntwf {c |}{col 14}{res}{space 2}  .810991{col 26}{space 2} .0389412{col 37}{space 1}   20.83{col 46}{space 3}0.000{col 54}{space 4} .7346677{col 67}{space 3} .8873144
{txt}W__ratio40~b {c |}{col 14}{res}{space 2}-.2912333{col 26}{space 2} .3486682{col 37}{space 1}   -0.84{col 46}{space 3}0.404{col 54}{space 4}-.9746104{col 67}{space 3} .3921438
{txt}W__ratio_5~b {c |}{col 14}{res}{space 2} 1.875044{col 26}{space 2} 1.011863{col 37}{space 1}    1.85{col 46}{space 3}0.064{col 54}{space 4}-.1081702{col 67}{space 3} 3.858258
{txt}B__ratio40~b {c |}{col 14}{res}{space 2} -1.83302{col 26}{space 2} .4333316{col 37}{space 1}   -4.23{col 46}{space 3}0.000{col 54}{space 4}-2.682335{col 67}{space 3} -.983706
{txt}B__ratio_5~b {c |}{col 14}{res}{space 2}-1.517599{col 26}{space 2} .8325484{col 37}{space 1}   -1.82{col 46}{space 3}0.068{col 54}{space 4}-3.149364{col 67}{space 3} .1141664
{txt}{space 7}_cons {c |}{col 14}{res}{space 2}-1.650608{col 26}{space 2} .9236549{col 37}{space 1}   -1.79{col 46}{space 3}0.074{col 54}{space 4}-3.460939{col 67}{space 3} .1597221
{txt}{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 4}/lnalpha {c |}{col 14}{res}{space 2}-3.297729{col 26}{space 2} .1813689{col 54}{space 4}-3.653206{col 67}{space 3}-2.942253
{txt}{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 1}{res}a_id        {col 14}{txt}{c |}
   var(_cons){c |}{col 14}{res}{space 2} .1816104{col 26}{space 2} .0323133{col 54}{space 4} .1281414{col 67}{space 3} .2573902
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}

{com}. *
. estat ic

Akaike's information criterion and Bayesian information criterion

{txt}{hline 13}{c TT}{hline 63}
       Model {c |}        Obs  ll(null)  ll(model)      df         AIC        BIC
{hline 13}{c +}{hline 63}
{ralign 12:{stata estimates replay model:model}}{col 14}{c |}{res}{col 16}       724{col 27}        .{col 38}-2352.639{col 49}    17{col 58} 4739.277{col 69} 4817.219
{txt}{hline 13}{c BT}{hline 63}
{p 15 21 2}
Note: N=Obs used in calculating BIC; see {helpb bic_note:[R] BIC note}.
{p_end}

{com}. 
. *
. *
. *
. *
. 
. 
. ** MODEL A3.4: RATIO OF OVER40 SUPERVISORS TO OVER40 NON-SUPERVISORS [RATIO OF ratio_40over_suplb TO ratio_40over_nonsuplb]: ONLY RANDOM INTERCEPT MODEL SPECIFICATION WITH BE & WE ESTIMATES FOR ALL COVARIATES [SANS YEAR UNIT EFFECTS & LN(TOTAL WORKFORCE)] --- FULL MODEL: ALL CONTROL COVARIATES ARE LNTEF & YEAR UNIT EFFECTS ***
. 
. xthybrid age_discrimination  ratio40suplb_nonsuplb    orgjustice_sem    nonprof40over_tr_lb     politicization_lb ratio_fsup_msup ratio_minsup_nonmsup lntwf ///
> yr2 yr3 yr4 yr5 yr6 yr7 yr8 yr9 yr10 ratio_55over_totallb if suffpaneldum!=0 , clusterid(a_id) vce(cluster a_id)  family(nbinomial) link(log) full ///
> use(ratio40suplb_nonsuplb   orgjustice_sem  nonprof40over_tr_lb    politicization_lb   ratio_fsup_msup ratio_minsup_nonmsup ratio_55over_totallb) 

{res}
{txt}{hline}
{p 0 8}Model {hi:model}{p_end}
{hline}

{txt}Mixed-effects GLM{col 49}{txt}Number of obs{col 67}={res}{col 69}       724
{txt}Family: {col 15}{res}negative binomial
{txt}Link: {col 29}{res}log
{txt}Overdispersion: {col 28}{res}mean
{txt}Group variable: {col 28}{res}a_id{col 49}{txt}Number of groups{col 67}={res}{col 69}        81

{col 49}{txt}Obs per group:
{col 63}{txt}min{col 67}={res}{col 69}         7
{col 63}{txt}avg{col 67}={res}{col 69}       8.9
{col 63}{txt}max{col 67}={res}{col 69}        10

{txt}Integration method: {col 21}{res}mvaghermite{col 49}{txt}Integration pts.{col 67}={col 78}{res}7

{col 49}{txt}Wald chi2({res}24{txt}){col 67}={res}{col 70}  1068.06
{txt}Log pseudolikelihood = {res}-2330.9667{col 49}{txt}Prob > chi2{col 67}={res}{col 73}0.0000
{txt}{ralign 78:(Std. Err. adjusted for {res:81} clusters in a_id)}
{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 14}{c |}{col 26}    Robust
{col 1}age_discri~n{col 14}{c |}      Coef.{col 26}   Std. Err.{col 38}      z{col 46}   P>|z|{col 54}     [95% Con{col 67}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 6}R__yr2 {c |}{col 14}{res}{space 2}-.0151638{col 26}{space 2} .0499128{col 37}{space 1}   -0.30{col 46}{space 3}0.761{col 54}{space 4}-.1129912{col 67}{space 3} .0826636
{txt}{space 6}R__yr3 {c |}{col 14}{res}{space 2}-.0563335{col 26}{space 2} .0559182{col 37}{space 1}   -1.01{col 46}{space 3}0.314{col 54}{space 4}-.1659312{col 67}{space 3} .0532641
{txt}{space 6}R__yr4 {c |}{col 14}{res}{space 2}-.1069798{col 26}{space 2} .0682401{col 37}{space 1}   -1.57{col 46}{space 3}0.117{col 54}{space 4} -.240728{col 67}{space 3} .0267685
{txt}{space 6}R__yr5 {c |}{col 14}{res}{space 2}-.1132099{col 26}{space 2} .0808123{col 37}{space 1}   -1.40{col 46}{space 3}0.161{col 54}{space 4}-.2715991{col 67}{space 3} .0451793
{txt}{space 6}R__yr6 {c |}{col 14}{res}{space 2}-.0078473{col 26}{space 2} .0765625{col 37}{space 1}   -0.10{col 46}{space 3}0.918{col 54}{space 4}-.1579071{col 67}{space 3} .1422125
{txt}{space 6}R__yr7 {c |}{col 14}{res}{space 2} -.008059{col 26}{space 2} .0788769{col 37}{space 1}   -0.10{col 46}{space 3}0.919{col 54}{space 4}-.1626549{col 67}{space 3} .1465369
{txt}{space 6}R__yr8 {c |}{col 14}{res}{space 2} -.110343{col 26}{space 2} .0811303{col 37}{space 1}   -1.36{col 46}{space 3}0.174{col 54}{space 4}-.2693553{col 67}{space 3} .0486694
{txt}{space 6}R__yr9 {c |}{col 14}{res}{space 2}-.0596718{col 26}{space 2} .0819652{col 37}{space 1}   -0.73{col 46}{space 3}0.467{col 54}{space 4}-.2203206{col 67}{space 3}  .100977
{txt}{space 5}R__yr10 {c |}{col 14}{res}{space 2}-.0751726{col 26}{space 2} .0944508{col 37}{space 1}   -0.80{col 46}{space 3}0.426{col 54}{space 4}-.2602928{col 67}{space 3} .1099476
{txt}{space 4}R__lntwf {c |}{col 14}{res}{space 2} .8034591{col 26}{space 2} .0363331{col 37}{space 1}   22.11{col 46}{space 3}0.000{col 54}{space 4} .7322474{col 67}{space 3} .8746707
{txt}W__ratio40~b {c |}{col 14}{res}{space 2}-.4563842{col 26}{space 2} .4418998{col 37}{space 1}   -1.03{col 46}{space 3}0.302{col 54}{space 4}-1.322492{col 67}{space 3} .4097235
{txt}W__orgjust~m {c |}{col 14}{res}{space 2} .1628627{col 26}{space 2} .2563088{col 37}{space 1}    0.64{col 46}{space 3}0.525{col 54}{space 4}-.3394932{col 67}{space 3} .6652187
{txt}W__nonprof~b {c |}{col 14}{res}{space 2}-.5007655{col 26}{space 2} 1.055961{col 37}{space 1}   -0.47{col 46}{space 3}0.635{col 54}{space 4} -2.57041{col 67}{space 3} 1.568879
{txt}W__politic~b {c |}{col 14}{res}{space 2}-1.965153{col 26}{space 2} 1.681054{col 37}{space 1}   -1.17{col 46}{space 3}0.242{col 54}{space 4}-5.259958{col 67}{space 3} 1.329653
{txt}W__ratio_f~p {c |}{col 14}{res}{space 2} .2714237{col 26}{space 2} .2457774{col 37}{space 1}    1.10{col 46}{space 3}0.269{col 54}{space 4}-.2102912{col 67}{space 3} .7531385
{txt}W__ratio_m~p {c |}{col 14}{res}{space 2}-.3552206{col 26}{space 2} .1108999{col 37}{space 1}   -3.20{col 46}{space 3}0.001{col 54}{space 4}-.5725804{col 67}{space 3}-.1378608
{txt}W__ratio_5~b {c |}{col 14}{res}{space 2} 1.983962{col 26}{space 2}  1.02224{col 37}{space 1}    1.94{col 46}{space 3}0.052{col 54}{space 4}-.0195915{col 67}{space 3} 3.987515
{txt}B__ratio40~b {c |}{col 14}{res}{space 2}-1.573944{col 26}{space 2} .3940768{col 37}{space 1}   -3.99{col 46}{space 3}0.000{col 54}{space 4} -2.34632{col 67}{space 3}-.8015677
{txt}B__orgjust~m {c |}{col 14}{res}{space 2}-1.580741{col 26}{space 2} .3449319{col 37}{space 1}   -4.58{col 46}{space 3}0.000{col 54}{space 4}-2.256795{col 67}{space 3}-.9046867
{txt}B__nonprof~b {c |}{col 14}{res}{space 2} .0256576{col 26}{space 2} .2943374{col 37}{space 1}    0.09{col 46}{space 3}0.931{col 54}{space 4} -.551233{col 67}{space 3} .6025483
{txt}B__politic~b {c |}{col 14}{res}{space 2} 1.831594{col 26}{space 2} .7617104{col 37}{space 1}    2.40{col 46}{space 3}0.016{col 54}{space 4} .3386693{col 67}{space 3} 3.324519
{txt}B__ratio_f~p {c |}{col 14}{res}{space 2} .2064888{col 26}{space 2}   .08966{col 37}{space 1}    2.30{col 46}{space 3}0.021{col 54}{space 4} .0307583{col 67}{space 3} .3822193
{txt}B__ratio_m~p {c |}{col 14}{res}{space 2}-.0502792{col 26}{space 2}  .051388{col 37}{space 1}   -0.98{col 46}{space 3}0.328{col 54}{space 4}-.1509977{col 67}{space 3} .0504393
{txt}B__ratio_5~b {c |}{col 14}{res}{space 2}-1.531423{col 26}{space 2} .6247404{col 37}{space 1}   -2.45{col 46}{space 3}0.014{col 54}{space 4}-2.755892{col 67}{space 3}-.3069543
{txt}{space 7}_cons {c |}{col 14}{res}{space 2}-1.984668{col 26}{space 2} .8924848{col 37}{space 1}   -2.22{col 46}{space 3}0.026{col 54}{space 4}-3.733907{col 67}{space 3}-.2354303
{txt}{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 4}/lnalpha {c |}{col 14}{res}{space 2}-3.338607{col 26}{space 2} .1748268{col 54}{space 4}-3.681262{col 67}{space 3}-2.995953
{txt}{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 1}{res}a_id        {col 14}{txt}{c |}
   var(_cons){c |}{col 14}{res}{space 2} .1229021{col 26}{space 2} .0273197{col 54}{space 4} .0794964{col 67}{space 3} .1900077
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}

{com}. *
. estat ic

Akaike's information criterion and Bayesian information criterion

{txt}{hline 13}{c TT}{hline 63}
       Model {c |}        Obs  ll(null)  ll(model)      df         AIC        BIC
{hline 13}{c +}{hline 63}
{ralign 12:{stata estimates replay model:model}}{col 14}{c |}{res}{col 16}       724{col 27}        .{col 38}-2330.967{col 49}    27{col 58} 4715.933{col 69} 4839.723
{txt}{hline 13}{c BT}{hline 63}
{p 15 21 2}
Note: N=Obs used in calculating BIC; see {helpb bic_note:[R] BIC note}.
{p_end}

{com}. 
. 
. 
. 
. log close
      {txt}name:  {res}<unnamed>
       {txt}log:  {res}C:\Users\jp18390\Dropbox\DISCRIMINATION PROJECT (U.S. FEDERAL AGENCIES)\Age Discrimination Project\PAR R&R\Statistics\Age_Discrimination.APPENDIX 3.ONLY LONG PANELS.06-02-2022.smcl
  {txt}log type:  {res}smcl
 {txt}closed on:  {res} 2 Jun 2022, 18:09:57
{txt}{.-}
{smcl}
{txt}{sf}{ul off}